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Abstract. We discuss basic concepts and properties of diffractive phenomena in soft 
hadron collisions and in deep-inelastic scattering at low Bjorken-x. The paper is not 
a review of the rapidly developing field but presents an attempt to show in simple 
terms the close inter-relationship between the dynamics of high-energy hadronic and 
deep-inelastic diffraction. Using the saturation model of Golec-Biernat and Wusthoff 
as an example, a simple explanation of geometrical scaling is presented. The relation 
between the QCD anomalous multiplicity dimension and the Pomeron intercept is 
discussed. 
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^ ■ 1. Introduction 

(N 

After nearly two decades outside of the mainstream of high-energy physics, the subject 
of diffractive scattering has made a spectacular come-back with the observation of Large 
Rapidity Gap events at the HERA ep collider and similar studies at the highest energy 
hadron colliders. It has become a field of intense research and many detailed aspects 
have been repeatedly reviewed JTJ,0]. 

To develop a phenomenology of very high energy scattering and diffraction, a field of 
research which originated in soft hadron hadron collisions, it is tempting and traditional 
to start from a i-channel formalism based on Regge theory. For deep-inelastic scattering 
(DIS) and hard diffractive processes this leads to the simple (and popular) Pomeron 
picture as first proposed by Ingelman and Schlein ||. Although Regge theory is perfectly 
valid and beautiful, based on very general properties of the scattering amplitudes, it is 
plagued by many problems in practical applications which, as happened in the past, 
severely limit its predictive power. Apart from the fundamental theoretical question 
how to derive the theory as a strong-coupling limit of QCD, the theory itself provides 
little insight into the relation between properties of the final states, the structure of the 
theory and the physical meaning of its parameters. This becomes particularly important 
when the formalism is used outside its traditional domain of application, such as in deep- 
inelastic scattering. 

In this paper, we shall mainly adopt an s-channel picture of diffraction. Diffractive 
scattering is explained by the differential absorption by the target of the large number 
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Figure 1. Kinematic variables (a) for the reaction e p — > e X; (b) for the semi- 
inclusive reaction e p ^ e N X with a rapidity gap. 

of states which coherently build up the initial-state hadron or (virtual) photon and 
scatter with different cross sections. Such an approach incorporates from the outset basic 
quantum mechanics and unitarity, and permits, at least conceptually, a unified treatment 
of hadron, and real and virtual photon scattering. It will allow us to appreciate the 
close inter-relation between the dynamics of high-energy hadronic and deep-inelastic 
diffraction (DDIS) at very small Bjorken-x and to understand that long-distance physics 
plays a very important role in both. It is also the approach used in the most successful 
of the present theoretical models JJ. 

The main thrust of the paper will be to argue that the physics can be understood 
on the basis of a surprisingly small number of dynamical ingredients. This in turn leads 
to a view of the collision dynamics which is simple enough to help develop intuition, 
provide physical insight and suggest fruitful avenues of research. The paper is mainly 
addressed to experimentalists entering the field and we hope it will broaden their view 
of the subject. No originality is claimed in the presentation of the material although 
any errors of interpretation should be attributed solely to the author. 

2. Preliminaries 

2.1. DIS kinematics and cross sections 

The standard kinematical variables to describe ep DIS are depicted in figure [l|a. The 
centre-of-mass energy squared of the ep system is s = (P + k) 2 , with P and k the initial- 
state four-momenta of the proton and electron (or positron), respectively. W, the CMS 
energy of the virtual-photon proton system, is given by W 2 = (P + q) 2 . The photon 
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virtuality Q 2 and the Bjorken variables x and y are defined as 



f = -Q 2 = (k — k')\ 
Neglecting the proton mass, one has 
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the latter expression being valid for x <C 1. 

For the "rapidity-gap" process presented in figure 03, and where a baryon with 
four-momentum P' is detected in the final state, one defines the additional variables 

Q 2 + M\ - t a Q 2 x 

~ = T 



t = (P-P 



p 



(3) 



Q 2 + W 2 ' ^ Q 2 + M 2 x -t 
The variable £ is the fractional energy-loss suffered by the incident proton. The variable 
(3 can naively be thought of as representing the fractional momentum carried by a struck 
parton in an object — Pomeron or Reggeon — carrying longitudinal momentum £, emitted 
by the proton and subsequently undergoing a hard scatter. For small \t\ one has 
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In strict analogy with the total ep cross section 
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which defines the structure function F 2 (x, Q 2 



is the QED coupling), the differential 



cross section for a semi-inclusive (SI) DIS process (figure |I]b) can be written as 
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Alternatively, in measurements of the diffractive (_D(3)) contribution to F 2 (x,Q 2 ) 
often uses the definition 
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dl3dQ 2 di /3 Q 4 [ ° 2(1 + R)_ 
replacing x by (3 in equation (|6|). R = a L /<j T is the ratio of the cross sections for 
longitudinally and transversely polarized virtual photons. Since y is usually small in 
experiment, R can be neglected. Equations @ and (^) are equivalent since they 
represent the same experimental data. From an experimental point of view, there is 
no a priori reason to prefer one over the other and both should be measured. For 
fixed (x, Q 2 ) (i.e. W fixed), the £ dependence of F 2 5/( ' 3 ' ) (£, x, Q 2 ) reflects that on M X - 
Alternatively, for fixed (f3,Q 2 ) (i.e. M x fixed) the W dependence of F 2 D ^(£, (3, Q 2 ) is 
explored by varying £. 

The structure function F 2 is related to the absorption cross section of a virtual 
photon by the proton, <t 7 * p . For diffractive scattering at high W (low x), we have 
similarly 

Q 2 
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2.2. Regge formalism 

2.2.1. Total and elastic cross sections. Since the Regge formalism is so often used in 
present analyses of diffractive HERA data, it is useful to recall here its main ingredients 
and predictions. For small-angle elastic scattering of two hadrons a and b at high s, 
dominated by Pomeron exchange, the Regge scattering amplitude (ignoring the small 
real part) takes the factorized form 

A$(s,t) = isP a (t)(s/s Y^)-^ b (t). (9) 

Here, s$ is an arbitrary mass scale, frequently chosen to be of the order of 1 GeV 2 . 
The dependence on the species of the incoming hadron is contained in the form factors, 
Pa,b(t), usually parameterized as an exponential cx exp(_B 0;a b t). a P {t) is the Pomeron 
trajectory. In its simplest version it is a Regge pole, with intercept a P (0) = 1 + e, 
slightly larger than 1; e controls the large- s or large W growth of the total and elastic 
cross sections. The name "Pomeron" was first used in ||, but first discussed by 
V.N. Gribov |J and later named after Y. Pomeranchuk J7|. The observed large- s- 
dependence of the cross sections can be accommodated with a trajectory ap(i) of the 
form 

a P (t) = a P (0) + a p t = 1 + e + a' p t. (10) 

In Regge theory a P (t), i.e. a P (0) and a' p , must be independent of the species of the 
particles colliding. We shall see that their meaning in terms of the particle production 
dynamics is, at least qualitatively, easy to understand. The energy dependences of the 
elastic and total cross section are given by 

do* /dt| t=0 = ^-iPa(0)(3 b (0)] 2 (s/s r, (11) 

of =f3 a (0)(3 b (0)(s/s y. (12) 

The meaning of a' p becomes clear if one considers da^/dt at small \t\, using an 
exponential approximation 

dotf/d* = — [(3 a (0)f3 b (0)} 2 e B{s)t (s/s o r = e B ^\ (13) 

l07T 107T 

with 

(14) 



B(s) = 2 (s 0; a + B . b + a'jp In . 
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The energy-independent terms B^ a ^ originate from the form-factors in equation (|9]). 
From pp data B . p ^2 — 3 GeV~ 2 . 

Equation (|HD shows that the forward elastic peak "shrinks" with energy: B(s) 
increases (here logarithmically) with s. In impact parameter (b) space (^j becomes 

A*(8, b) = i ^(^(0) (^0) e e -g»/2B W (15) 

The transverse size of the "interaction region" is Gaussian with B(s) = (b 2 ). 

The collision can be visualized in the impact-parameter plane, figure |2|. The 
scattering profile is a disc with a 6-dependent opacity, the mean radius of the disc 
being proportional to ^jB{s). B(s) contains an s-independent and a lns/so term. 
The radius expands with s with a rate of growth determined by a' p , estimated to be 
~ 0.25 GeV -2 [U (for a much earlier but similar estimate see ||). At the same time, 
the opacity at fixed b is likely to increase too. In a perturbative QCD picture, this 
corresponds to an increase of the gluon density in the target or, in the proton rest frame, 
with an increase of the interaction probability (symbolized by increasing blackness in 
the figure). Since the latter cannot exceed unity, it follows that also the gluon density 
cannot rise indefinitely. 

Concerning the meaning of a F (0), it will become clear in section |6.2| that the 
increase of <tt with s can be attributed to an increase of the number of (wee) partons 
in projectile and target. This will allow us to relate a F (0) to the wee-parton density in 
rapidity or, more generally, to the QCD multiplicity anomalous dimension. 



2.2.2. Triple-Regge parameterization of the reaction a + b — > c + X . Regge theory can 
be generalized to inclusive reactions. The invariant cross section of an inclusive process, 
a + b — > c + X, can be expressed in terms of the Regge- Mueller expansion which is based 



on Mueller's generalized optical theorem [[10]. This states that an inclusive reaction 



ab — > cX is connected to the elastic ("forward") three-body amplitude A(abc — > abc) 
via 

E ^-^-(ab — > cX) ~ - DiscM 2 A(abc — > abc), (16) 
dp 6 s 

where the discontinuity is taken across the M\ cut of the elastic Reggeized amplitude. 
For the triple-Pomeron diagram, valid in the (diffractive) region of phase space where 
the momentum fraction of particle c is near one, or s > W 2 ^> (Mj^Q 2 ) 3> (|t|,m 2 ), 
equation (|i~6|) has the approximate form 

„ d 3 cr 1 d 2 cr s d 2 cr „ s 

E ^ = -,^-,wm = mt) - aM ' (17) 

where t is the four-momentum transfer squared. The "flux factor" f(£,t) is given by 
f(£,t) = N F 2 (t) fl 1 " 2 ^*)]. (18) 
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In the above equations, N is a normalization factor, F(t) the form-factor of the ppIP- 
vertex; o~$> p can be interpreted as the Pomeron-proton total cross section. Assuming 
0>p ~ (M x ) e , and using equation (|I0|), (|T7|) and (|T3D imply that 



rfV s £ (M 2 



^_ fig) 

£(l+e) £(l+2e)" ^ 



The two expressions on the right-hand side equation ( |T9| ) are equivalent. However, they 
show that the model predicts slightly different £ dependences if either s is kept fixed and 
Mx varied, or if M x is fixed and s varied. Nevertheless, since e is small, both expressions 
in equation flT9| ) show that the triple-Regge IPIPIP contribution in the region £ <C 1 is 
of the generic form 



d 2 a 



dt d£ 



with S < 1 (20) 



t=o 



£(1+5) 

and predicts a "universal" l/£ dependence as long as 5 is universal. Regge theory also 
implies that both a e \/ai and cx d /<7t increase as s t . This eventually leads to violation of 
unitarity since e is found to be positive. The total diffractive cross section, cr^, grows 



as s 2t . 



Although the applicability of Mueller's optical theorem to reactions with (far) off- 
shell particles has not been proven, it is very frequently used as the starting point 
in analyses of diffractive phenomena in 7*p scattering. In that case, s in the above 
equations has to be replaced by W 2 or 1/x if Q 2 is fixed. 

2.2.3. Problems. In spite of the elegance of the Regge approach, it has been known 
for a long time [|ll|] that the theory with a "super-critical Pomeron": ajp(0) = 1 + e 



(e > 0), is plagued by unitarity problems as s — ► oo which are especially severe for 
inelastic diffraction: i) the power-law dependence, <tt oc s e violates the Froissart-Martin 
bound [[TIJ; ii) the ratio u c \/a^ oc s e /lns eventually exceeds the black-disk geometrical 
bound [a e i < \&t)\ hi) the ratio (t°/<tt increases as s e . This disagrees with experiment 



not only for hadron collisions [13|, but also for deep- inelastic diffraction, where the ratio 



Or/or is found to be essentially independent of W |14|,|l5|l (see section |4.3|) . 



3. Experimental results on total and elastic cross sections 

3.1. Energy dependence of hadronic total cross sections 

The s-dependence of total hadron-hadron cross sections, <tt, has been measured for many 
combinations of hadrons. Above ~ 20 GeV all hadronic cross sections rise with s. This 



was first discovered for K + p collisions in 1970 at the Serpukhov accelerator |fL7|1 . The 



rise of the pp total cross section was first observed at the ISR |L8| and later confirmed 



at Fermilab |1| . A compilation of pp, pp and ^p data is shown in figure The solid 
lines are fits which include a component decreasing rapidly with s and a second rising 
component which persists at high energies. In |16|] it was observed that all measured 




Figure 3. Total cross sections measured in hadronic scattering as a function of the 
centre-of-mass energy for pp, pp, ir^p scattering. The cross sections show a "universal" 
rise at high energies of the form a ~ s 08 |l6| ]. 



hadron-hadron (and jp) cross sections grow in an similar way at high s. An economical 
parameterization is a sum of two power-law terms in s 

a T = Xs e + Y s e \ (21) 

where the constants X and Y depend on the reaction. This obviously is inspired 
by Regge theory, the two terms in equation (|2T|) corresponding to Pomeron and 
"normal" Regge (Reggeon) exchanges, respectively. The value of e is not very precisely 
established. Various recent "global" fits find the data to be compatible with e in 



the range 0.08 - 0.1 pi, EDI, pTJ ; e' is found to be 0.45 pL Global fits to total 



elastic and diffractive cross sections performed much earlier yielded similar values for 



One should also note that the present data cannot discriminate between 
"simple-pole" fits inspired by a Regge-model of t-channel exchanges leading to a power- 
law dependence, and equally valid fits to log 2 s and log s (or, for that matter, e^ logs ) 
functional dependences. 

Although the significance of equation ( EI|) has been over-emphasized, the 
"universal" high-energy behaviour of the total hadronic cross sections is an important 
observation which calls for deeper understanding. It also raises the question (not 
addressed in Regge theory) which particular final states are responsible for this increase. 

As already mentioned, the single-Pomeron exchange amplitude violates unitarity 
thus indicating an inconsistency of this model. The simplest way to overcome this 
problem is to introduce multiple Pomeron exchanges (or multiple interactions) in a 
single scattering process, as shown in figure [| p4| . The total amplitude can then be 



written as the sum of n-pomeron exchange amplitudes A^ n \s,t). For each n-pomeron 
graph one can define a theoretical "total" cross section applying the optical theorem to 
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-. oo 
a (n) = (_l)n+l£ Imm ( A (n)) ? ^ = ^(_1)™+V (n) . (22) 

S n=l 

As a simplified model consider only the first two graphs shown in Fig. [|, assuming 
a ( n ) <^ 4cr( 2 ) < orW with n > 2. Then, the total cross section becomes a to t = 0"^ — 0"^ 2 \ 




(a) (b) (c) 

Figure 4. Hadron-hadron scattering via pomeron exchange: (a) one-Pomeron, (b) 
two-Pomeron, and (c) three- Pomeron exchange graphs. From J&jj/ . 



where and are the cross sections of the one- and two-pomeron exchange graphs, 
respectively. The energy-dependence of the two-pomeron cross section is directly linked 
to that of s £ and turns out to be The two-pomeron cross section 

grows faster with energy than the one-Pomeron cross section. Since its contribution is 
negative, this leads to a weaker energy-dependence of the total cross section than in 
the single-Pomeron exchange model and a smaller effective Pomeron intercept. It also 
breaks Regge factorization. 

Interestingly, according to the Abramovski, Gribov, Kancheli (AGK) cancellation 
theorem, the contribution of the two-pomeron graph to the inclusive inelastic 



single-particle cross section vanishes ||26|| . Analogously, the factorization violating 
contributions due to multi-Pomeron exchange graphs cancel out in all orders. This 
means that only the one-Pomeron graph determines the inclusive particle cross section 
in the central region. Consequently, a study of the energy-dependence of the single- 
particle inclusive spectrum should allow to mesure the value of the Pomeron intercept in 
soft hadronic interactions, in a way which is unaffected by multi-Pomeron (or screening) 
effects. 

The results of such an analysis is shown in figure |] f25| . The authors use a double- 
Regge expansion, valid at high energies and in the central region of centre-of-mass 
rapidity (y — 0), which predicts the energy-dependence 

^ = a PP s A + a RP s^-V + a RR (23) 

dy y=Q 

where the a-parameters are Reggeon couplings and 1 + A is the value of the Pomeron 
intercept, unaffected by multi-Pomeron absorptive effects. The fit yields A = 0.170 ± 
0.008, for negative particle (c~) production and A = 0.167 ± 0.024 for K° s inclusive 
production. As anticipated above, this is substantially larger than the effective intercept 
~ 0.08 deduced from the s-dependence of hadron-hadron total cross sections which is, 
however, affected by the rescattering contributions. 
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Figure 5. Cross sections of negatively charged particles and Kg in the reactions 
pp — > c~ + X, pp — > c~ + X , pp — > Kg + X and pp — > Kg + X in the central centre-of- 
mass rapidity region, (da/dy) y —o. c~ stands for a negatively charged hadron. The solid 
curves are fits with the double-Regge expression ( p3| ) with a super-critical Pomeron: 
as A +bs( 2A ~ 1 ^ 4 +cs~ 1 / 2 yielding A ~ 0.17 for all reactions. The dashed lines represent 
the s A term. For references to data see pal. 



In |27] it is argued that the "bare" value of A is still larger, since renormalization 
effects induced by Pomeron-Pomeron interactions lower its effective value. The 
correction is estimated to be ~ 0.14. In all, this implies that the bare Pomeron intercept 
could be as large as 1.3 and thus comparable (see below) to what is measured in deep- 
inelastic scattering. For the latter process, absorption effects due to multi-Pomeron 
exchange are expected to be much smaller than in soft hadronic collisions, due to the 
short interaction time, and to diminish with increasing Q 2 with the result that the 
Pomeron intercept measured in DIS could come close to that of the bare Pomeron 
"active" in soft hadron collisions. 

Let us note also that the parameterization (E3[) predicts cross sections for negatively 



charged particles and K$ of 251 ± 26 mb and 25 ± 7 mb, respectively, at the LHC. 
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Figure 6. 7*p cross section as a function of W 2 at various Q 2 , shown on the left 
side together with the scale factor applied to the data for better visibility. The full 
lines show a QCD-fit p8| , the dashed lines are a fit by the Golec-Biernat Wusthoff 
saturation model pal . From [BOl . 



3.2. The 7*p total cross section at HERA 

The measurement of the total 7*p cross section as a function of Q 2 and W is one of the 
major achievements of the experiments at HERA. Some results are shown in Fig. |6] ||32j| . 
Remembering that 1 / \/~Q 2 ~ = i? 7 * determines the transverse distance which the photon 
can resolve, we note that for small Q 2 (large FLy*) the cross section has a hadron-like 
increase with W: the photon acts like a hadron. With increasing Q 2 , the rise with W 
becomes stronger: the photon shrinks and becomes more and more point-like. 

Parameterizing the jy-dependence as a] [ ~ (W 2 ) Xtot , one obtains the results shown 
in figure Within the measured range, Xtot increases linearly with log (Q 2 ) from a value 
~ 0.08 at low Q 2 , the same as in hadron-hadron interactions, to ~ 0.35 at the highest 
Q 2 . These data were also analysed in |[33|| . If interpreted in terms of Regge exchanges, 
it is clear that for 7*p collisions, "universality" of the trajectory parameters no longer 
holds: ajp(O) depends on Q 2 , and a continuous transition is seen between the soft regime 
and that where a "small-size" 7* hits a proton. The dynamics evolves in a continuous 
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Figure 7. The exponent A tot in a] [ ~ (t¥ 2 ) Atot , versus Q 2 . The full line 
shows a fit to the form indicated with a — 0.0481 ± 0.0013(stat) ± 0.0037(syst) and 
A = 292 ± 20(stat) ± 51(syst) MeV. Extrapolation to Q 2 = 0.48 GeV 2 gives a value of 
0.08l3l. 
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manner. Evidently, as is well-known, the results for "small-size" virtual photons can be 
(partly) interpreted in terms of perturbative QCD radiation and the familiar parton- 
density evolution equations. 

3.3. Elastic scattering and forward slope 

3.3.1. Hadron hadron interactions. Figure |8|a shows data on the forward elastic slope 
in pp and pp interactions. The shrinkage of the diffractive peak with -y/i, expected from 
Regge theory is clearly seen. Expressed in geometrical or optical terms, the "effective 
interaction radius" of the proton becomes larger with increasing s, as schematically 
illustrated in figure |2|. 

The values of the slopes are in rough agreement with what is expected for (optical) 
diffraction on a "black" fully absorbing disk of radius R for which B = R 2 /A. For 
a proton with R w l/wv (m n is the pion mass), B is expected to have a value of 
13 GeV -2 which compares well with the data. However, for scattering on a black disk, 
Oci/°T = 1/2, whereas experiment, figure ||b, shows a value between 1/5 and 1/4 at high 
s. This means that the proton is semi-transparant, even at zero impact parameter as 



shown experimentally in pl |. Indeed, since the wavefunctions of the hadrons entering 
the collision are a superposition of states, some will be fully absorbed, while others will 
pass through almost unaffected. This agrees with the idea of color transparency in QCD 



see section |6.3.2| ). Such a mixture of states with very different absorption probabilities 



will be essential for inelastic diffraction to occur, see section p7| 
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Figure 8. (a) slope parameter B(s), (b) ratio of elastic to total cross section versus 
a/s for pp/pp interactions. The solid lines are Regge fits. For details see pl|. 
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Figure 9. W dependence of the cross section a{^*p — ► p°p) for various Q 2 values. 
The data for Q 2 < 1 GeV 2 obtained previously |35| are also shown. The solid lines 
show a fit with <J 7 * p ^ p o p ~ W s . The shaded area indicates normalization uncertainties 
due to proton dissociation background. From [f36|. 



3.3.2. Real and virtual photon quasi-elastic scattering. Among the many results now- 
available (for a review see fl37|), figure |9] shows, as an example, DIS measurements of 
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the W- dependence of elastic p° electroproduction as a function of Q 2 . For each Q 2 
interval, the cross section is assumed to be of the form W s . In the same manner as 
for the 7*p total cross section, the data suggest a marked increase of 5 when Q 2 enters 
a regime where pQCD becomes relevant. However, the errors remain sizable and, in 
PU-regions where the DIS data overlap (3.5 < Q 2 < 13.0 GeV 2 ), the Q 2 -dependence 
of S is statistically not yet very significant. Measurements at larger Q 2 are needed to 
clarify this important issue. 
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As to the shape of the diffractive peak, figure [T^ shows a compilation |30| of the 

2 



slope B, at fixed PV, as a function of an effective scale, Q 2 



Q 2 + my, for various 



vector mesons with mass my. 



The slope becomes smaller with increasing Q 2 



In 



the photoproduction region, Qjjf = 0, the slopes for p and u> are quite similar to those 
observed in proton-proton scattering, see figure §. At higher Qlff, they are considerably 
smaller, approximately half of that observed in proton-proton scattering. The effective 
interaction region reduces to about that of a single proton, as expected for a projectile 
which becomes more point-like as Q 2 grows. At the same time the total cross section 
itself grows faster with W than in hadron interactions. 

For elastic J/ip photoproduction, ZEUS recently measured (see P5[ ) the differential 
cross-section, da/dt oc (W) 2oiF ^~ 2 , in the energy-range 20 < W < 290 GeV, 
and \t\ < 1.25 GeV 2 . From the t-slope in bins of W, yielding e.g. B = 
4.3 ± 0.08 (stat)t°;^ (syst) GeV" 2 at W = 90 GeV, one derives a P (0) = 1.201 ± 
0.013 (stat)l^?? (syst) and a' p = 0.126 ± 0.029 (stat)l^ (syst) GeV" 2 . The latter 
value implies that shrinkage is smaller than in soft hadronic collision but not negligible. 
This was predicted in |3(| . 
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A recent first-time ZEUS measurement of the leading trajectory parameters 
from exclusive p production in DIS, 7*p — > p°p, (1 < Q 2 < 40 GeV 2 ), yielded 
a P (0) = 1.14 ± 0.01 (stat)lg;^ (syst), = 0.04 ± 0.07 (stat)^ (syst) GeV -2 . While 
not conclusive, given the errors, this measurement also suggests a smaller value of a' p 
than that of the "soft" Pomeron (a' F « 0.25 GeV -2 ). 

From a measurement of the spin-density matrix of the p° decay, ZEUS also 
extracted o"l/ct^, the ratio of the cross section for longitudinally and transversely 
polarized 7*, as a function of Q 2 and W. The ratio strongly increases with Q 2 but 
is found to be independent of W, a somewhat surprising result given the expectation 
that at large Q 2 the average transverse size of the longitudinally polarized 7* is much 



smaller than that of a transversely polarized 7* [40 



3.3.3. Brief summary. Although more precise measurements are evidently needed, and 
forthcoming, the present data on total and elastic differential cross sections suggest 
a clear trend. As for real hadrons, for near-on-shell photons fluctuating into light 
vector mesons, (p, <f>) and which have large (order 1 fm) transverse extensions, inversely 
proportional to the Compton wavelength of the light quarks in the meson, the effective 
Pomeron trajectory a P (t) is close to that of soft collisions. For heavier vector mesons 
(e.g. J/V), which are characterized by a smaller transverse size, or in DIS, present data 
provide some indication for a weaker shrinkage, with a' p smaller than the "soft" value 
0.25 GeV -2 . At the same time, the effective intercept a P (0) grows with decreasing size 
i? 7 *. The transition from the soft hadron-like regime to DIS is a smooth one. 



4. Inelastic diffraction 

4-1. Experimental signatures 

In contrast to forward elastic scattering, which beautifully reflects the wave-nature of the 



particles, the phenomenon of diffraction dissociation, predicted by Good and Walker [II 
has no classical analogue. For hadron-hadron scattering, it corresponds to quasi-elastic 
scattering between the two hadrons, where, in single diffraction, one of them is excited 
into a higher mass state retaining its quantum numbers. This coherent excitation, 



illustrated in figure |TT| for single-diffraction, requires not only small transverse (APr) 



but also small longitudinal (AP^) momentum transfer. This leads to the coherence 



condition (see e.g. p2|,|43 



^M<I^~o.l-0.2. (24) 



s m p 



The coherence condition arises from the need to conserve the coherence of the quasi- 
elastically scattered target and implies that the diffractive mass Mx cannot be too large. 
For zero-angle production the minimum four-momentum transfer at which the mass Mx 
can be produced is \t min \ = [(M x — m 2 )/2p] 2 , with p the incident momentum in the 
target rest frame. In the transition, the wavenumber k of the incident hadron varies 
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duced hadrons, M, is denoted 

by M x in the main text. Figure 12. Topology for pp -»■ pX Q. 

From Q. 



by an amount Ak oc ^/|t m m|. The condition of coherence follows from the requirement 
that the wavenumber changes little during the passage through the target, so that the 
waves describing the target before and after the interaction can stay in phase. For 
DIS kinematics, the minimum value of t required to produce a given Mx from a target 
with mass itit is |t m m| — rn\[M\ + Q 2 ) 2 /W A . For a typical hadronic radius of 1 fm, 
M\ < 0.2 W 2 . 

The generic topology of a single-diffractive (here pp) event is illustrated in figure |T^ . 
The upper-limit on Mx implies that the diffractive hadronic final states exhibit a large 
rapidity gap between the quasi-elastically scattered proton and the dissociation products 
X of the p. The width of the gap in (pseudo-)rapidity space measured from the rapidity 
of the initial-state proton is At] ph In |. In collider experiments diffractive events are thus 
identified either by detecting directly a "fast" ( "leading" ) proton in a spectrometer, by 
the presence in the main detector of a large rapidity region devoid of hadrons (a rapidity 
gap), or by exploiting the characteristic l/M x dependence of diffraction. 

Naively, the interaction is often viewed as proceeding via the emission from 
the proton of a Pomeron, a colorless object with vacuum quantum numbers which 
subsequently interacts with the p. In QCD such an object, if it were to exist as a 
physical entity, must be a colour-singlet composed of quarks, antiquarks and gluons. It 
will become clear later, however, that such a picture is an unnecessary and probably 
misleading simplification of mechanism behind diffractive physics. 



4-2. Hadron hadron inelastic diffraction 

Evidence for an important diffractive component in the inclusive reaction p + p — > 
p + X, with excitation of large masses, was first established at the ISR by the CHLM 
collaboration j|6 . Figure [TJ| shows single diffractive pp cross sections from low to high 
s. The diffractive enhancement becomes less and less prominent as s decreases, in line 
with the previous discussion about the need to maintain coherence of the target. The 
Mx-spectrum drops rapidly in the resonance region. Beyond that it levels-off and shows 
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Figure 13. Single diffractive pp cross sections p3|). The figure shows how the 
characteristic \jM\ (Regge) behavior of diffraction becomes manifest as s increases. 
The arrows indicate the low mass (Mx ~ 2 GeV) resonance region. From . 
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an approximate 1/M X dependence. 

A compilation of measurements now plotted against M x , is shown in figure [TJ 



for and single diffractive cross sections at t — —0.05 GeV 2 (for earlier compilations 
see The distribution falls as 1/ (M x ) 1+A over the entire M x region. Quite 

remarkably, it is independent of s over five orders of magnitude. The data are consistent 
with the same value of a P (0) — 1 = e = 0.104 (denoted A in the figure) as that extracted 
from the fit in [^l|] to total and elastic cross sections data. 

The 1/M X scaling shown by the data in figure [14] implies that, since e is small, the 
rapidity-gap distribution, dV^/dt dArj is nearly independent of s. In Regge models of 
diffraction, this distribution is related to the "Pomeron flux", see equation ([i~8T). Such a 
weak energy-dependence must reflect a fundamental, but not yet understood, property of 
the baryon "re-formation" process in the final state. A strict energy-independence would 



be consistent with short-range order |48|, a basic property of multiparticle production. 



4-3. Inclusive diffraction at HERA 

In DIS at small x measured at HERA, inelastic diffraction occurs at a rate of 0(10%) 
of all events ||49 |. Although it surprised many in the pQCD community, it had been 



anticipated even before the advent of QCD [|5(J. It was also predicted from Regge 



theory |5T| . The occurrence of such diffractive events, also called "Large Rapidity Gap 



events" are indeed difficult to understand in the parton picture on the basis of pQCD 
alone. 

The experimental effort at HERA has concentrated on measurements of the 
diffractive part, F^ 3 \ of the structure function F 2 , equation (^). The data have been 
reviewed on many occasions and details can be found in [1], |52|.[53| . 



New preliminary HI 1997 inclusive diffractive data [|Tj| have been used to extract 
a P (0) from the ^-dependence of (/3, Q 2 , x P ) with much increased precision, 

yielding a P (0) = 1.173±0.018(stat.)±0.017(syst.)^ Q 3 5(model). This value is not much 
higher than a P (0) ~ 1.1 in soft processes. As figure [15] shows, there is no evidence for a 
systematic variation with Q 2 . The data further suggest that the effective intercept for 
°t(7*p) is larger than that of the diffractive contribution at high Q 2 . 



Another striking HERA result, first obtained by ZEUS []TJ], is illustrated in figure [16] 
and in figure |T7| which shows more recent HI measurements |T3 . For Q 2 and M x (and 



thus (3) fixed, the relative rate of diffractive events is nearly PU-independent, except at 
very small (3. Standard triple-Regge theory, without multi-Pomeron exchange, predicts 
an increase as (W 2 ) e , in clear disagreement with data. 

As discussed in f55|, no adequate explanation within purely pQCD of the constancy 
of the mentioned ratio is known at present. The authors conclude that the non- 
perturbative QCD contribution to diffractive production is essential. Indeed, constancy 
of the ratio is obtained quite naturally in the quasi-classical gluon field approach (see 
Buchmuller in 0). It is also correctly predicted in the GBW- model ||29|| . There it is a 
consequence of the basic assumption that the cross section of the system radiated off 
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Figure 15. HI: effective value of a F (0) as function of Q 2 : • extracted from F 2 
■ from a fit of F 2 {x, Q 2 ) to the form ar A (Q 3 ); k HI 1994 data @. From p[. 
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Figure 16. Ratio of diffractive and total cross sections at fixed values of Q 2 , for 
different regions of the invariant diffractive mass Mx- The lines are predictions of the 
saturation model [fEl . 



the 7* partonic fluctuations saturates once this system has acquired a large transverse 
extension and is thus non-perturbative. 

4-3.1. An interlude: Leading baryons, the energy-loss spectrum. Inclusive diffraction 
shows a fractional proton energy-loss spectrum of the form given in equation For 
e = this intruigingly resembles a soft bremsstrahlung spectrum. Early QCD-models 
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Figure 17. Recent HI measurements of p D< - 3 \ the ratio of the diffractive to the 
inclusive cross section versus W. Data at (3 = 0.9 are scaled by a factor of five [15j. 



for diffraction have been proposed based on this analogy which indeed predict a \jM\ 
spectrum [fj£|. In this context it is interesting to recall a calculation of the leading- 
particle (LP) energy loss using a QED soft radiation analogy by Stodolsky [j5{| . 



Assume that a "leading" particle loses energy analogously to an electron which emits 
soft photons via bremsstrahlung. If ( is the total energy lost by the incident hadron 
which has initial energy Eq, the probability to radiate N particles of total energy £, 
with Ni of them having energy u>i, is given by 

M0 = J2 P ( N ^ P ( N ^ ---HC-NxUt- N 2 uo 2 -■■■) 

Ni,N 2 ,... 

x 5(N - Ni - N 2 ), (25) 

where 

f w- J exp (- W (26) 

is the Poisson probability (valid in QED for soft radiation) for Ni emissions in the energy- 
interval u>i, u>i + did, with [(dN / cko)did} their mean number. Setting dN/dcu = \/u and 
summing over all N one finds after a lengthy calculation the surprisingly simple result 

f(z) = ~ = X(l-\z\) x - 1 ] (27) 

with z the fractional energy E / 'Eq of the proton. For A 1 one obtains a flat 
distribution. 
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Figure 18. Inclusive spectrum for the reactions pp —> pX and pp — > pX (insert). 
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Figure 19. The fraction of leading protons, measured with the ZEUS Leading Proton 
Spectrometer, as a function of Xl = 1 — £ in bins of a; and Q 2 . From [|57f . 
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The previous calculation assumed Poisson emission which disagrees with 
experimental observations (and pQCD predictions). Generalizing to an emission process 
where the multiplicity distribution obeys Koba-Nielsen-Olesen (KNO) scaling (valid in 
pQCD (N)P(N) = V(N/(N)) with V(u) energy-independent, one finds §2| 



1 r\rr /"°° 

- — = J d(A/A)*(A/A) A(l - z) 1 -". (28) 

The mean number of radiated objects is (N) = A ln(s/so) with so a scale parameter. 
A flat z spectrum is recovered for A ~ 1. The parameter A is the mean number of 
emitted objects per unit of rapidity. For comparison with experimental data one has to 
assume that these objects (they were called "clusters" in ancient times) are resonances 
or higher-mass states decaying on average into two or three final-state particles. A 
density A ~ 1 is therefore a reasonable number. 

In the limit z — > 1, and for ty(u) ~ u 13 near u — 0, one obtains 

^~ r <' 3+2) r^ H i-Vr" - (29) 

Thus, besides being flat away from z = 1, the spectrum develops a diffractive- like peak 
at large z. Ignoring the logarithmic factor, this result coincides with equation (pC|) 
for ajp(0) = 1. In fact, in triple-Regge language, the full expression, equation (|29|), 
corresponds to a Pomeron cut, and not a simple Pomeron pole, in line with general 
theoretical expectation. In this model, the enhancement near z = 1 is due to low- 
multiplicity events. The detailed shape of the spectrum is, therefore, determined by 
that of the KNO function at small values of u. 

We find it remarkable that the quite simple and reasonable assumptions leading 
to equation (EH) are sufficient to capture essential aspects of the LP spectrum and its 
'diffractive limit", z — ► 1. If both A and ^ vary slowly with energy, the same will 
hold for the LP spectrum and for the diffractive peak and is not in disagreement with 
experiment. 

The "flatness" of the leading proton spectrum is well-known from hadron-hadron 
collisions. An example for pp interactions is shown in figure |18|. The same flatness is 
seen in DIS data (an example from ZEUS is shown in figure [19|). Although the spectra 
are found to be independent of x and Q 2 in the DIS regime, a small but significant 
increase of the rate with Q 2 is now seen in the low-Q 2 region | 57fl . 

Whereas diffractive data at very small £ and so-called Leading-Baryon data 
at larger £, outside of the diffractive region, are usually analyzed separately, our 
previous discussion argues in favour of combined analyses of such data. This has 
recently been done by the authors of who combined diffractive structure function 
measurements with Leading-Proton and Leading- Neutron results from HI. Results 
are shown in figure |20] |15| which displays £ F 2 LP and £ as functions of £. A 
combined triple-Regge fit, including Pomeron, Reggeon and pion exchange contributions 
yields a P (0) = 1.250 ± 0.023 and a Reggeon trajectory compatible with / 2 exchange: 
«m(0) = 0.770 ± 0.030. Note the somewhat larger value of 0:^(0) than the recent HI 



measurement quoted in section [O 
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Figure 20. Diffractive (Diff, open circles and open squares) and leading proton (LP, 
black triangles) structure function data vs. £, for fixed x and Q 2 . The figure combines, 
in each plot, the diffractive and leading proton HI data with similar values of x and 
Q 2 . The quoted Q 2 values are those of the diffractive data; the corresponding Leading 
Proton values are Q 2 = 4.4, 7.5, 13.3 and 28.6 GeV 2 . The black circles and black 
squares are data with Mx < 2 GeV. The plotted curves represent a global fit: total 
(solid line), pomeron (dotted line) and reggeon plus pion contributions (dashed line). 
For details ands references see lEal. 



5. Unitarity 

The importance of unitarity is not always sufficiently appreciated. This section is 
therefore devoted to a description of its main aspects and implications. 

The unitarity of the scattering matrix, T, implies close relationships between total 
cross sections, the elastic scattering amplitude and the amplitudes of inelastic final 
states. The unitarity relation between states \i> and |/> reads 

21m <f\T\i>= <f\T + \e><e\T\i> + < f\T + \n><n\T\i> } (30) 

|e><e| \n><n\ 

where Yl\ e ><e\ stands for summation and integration over all possible elastic 
intermediate states |e>. The second term is the contribution from all possible inelastic 
states; \i> is the initial and |/> an arbitrary final state. 

For forward elastic scattering, t = (\i >= \f >) equation ( |3~0"|) immediately leads to 
the optical theorem. However, the relation has much wider consequences since the state 
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| / > can be any state. It shows that the imaginary part of the amplitude of any particular 
final state </|T|«> in general receives contributions from all other final states. Such 
"unit arizat ion effects" will be small only if the "overlap" (< f\T\n >< n\T\i >) of the 
states |i>, \ f> with the states |n> happens to be small. This will, therefore, depend 
crucially on the topology in momentum space of the inelastic states and on the phases 
of the amplitudes. 

The two terms on the right-hand side of equation ( |30l) are called the elastic and 



inelastic overlap functions, respectively, and were first introduced by Van Hove |63 
For elastic scattering (and neglecting the real part), we see that the amplitude can, in 
principle, be calculated from the knowledge of the inelastic final states. This is the 
so-called s-channel approach to diffractive scattering. It provides an alternative to the 
t-channel approach in which the diffractive amplitudes are analyzed in terms of their 
singularities, poles and cuts, in the complex angular momentum plane. 

An important result is obtained (valid only at high s) when equation (|30D is written 
in impact-parameter (6) space|]]. Using angular momentum conservation one finds 

21m Ad{s,b) = \A el (s,b)\ 2 + G in (s,b). (31) 

Here A e i(s,b) is the elastic amplitude; G in (s,b), the inelastic overlap function, is the 
contribution from all inelastic channels. From equation ( |3l"l) follows that Im A e i(s,b) 
at impact parameter b is generated by the absorption into the inelastic channels at the 
same impact parameter: "unitarity is diagonal in 6-space" . 

For ReAei = 0, equation (|3l"D can be solved easily for A e i if Gi n (s,b) is known. 
Alternatively, knowledge of A e i(s, b, t) can be used to determine G in (s, b) (see e.g. p4f ). 
For DIS, it is presently unknown but of great interest. 

Equation (|3T| ) has the general solution 

G m (s,b) = l-e- n ^ (32) 

A d (s,b)=i {l - e-^^ b \ } . (33) 

The "opacity" function or eikonal, £l(b), and the phase x(s,6) are arbitrary real 
functions. The former has a simple meaning: exp [— fl(b)} is the probability that no 
inelastic interactions with the target occur. We further have the general relations 

cTei{s) = J d 2 6|^(s,6)i 2 , (34) 
a T ( s ) = 2 J d 2 b Im A e i(s,b), (35) 
<r in {s) = J d 2 b [21m A el (s,b) - \A el \ 2 (s,b)] . (36) 

5.1. Elastic diffraction and shrinkage 

For scattering on a proton, absorption into inelastic channels will be most important 



for values of b smaller than the proton radius. From equation (|3~T|) follows that this 



|| For a recent mathematical discussion of the validity of this transformation at finite energies, and 
further references, see |p3|. 
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will generate a large imaginary elastic amplitude at the same impact parameter. The 
impact parameter profile will be maximum at b = 0, where absorption is strongest, 
and the elastic differential cross section, da e \/dt, sharply peaked at t = 0, its width 
reflecting transverse extension of the effective interaction region. The experimental fact 
that Re A e \ is small at high s implies that elastic scattering can indeed be considered as 
the "shadow" of the inelastic channels. 

The physical meaning of the slope B(s) can also be understood from the shape of 
Gi n (s, b) and equation (|30D . Indeed, Gi n (s, b) is a measure of the overlap of the amplitude 
of a given final state with the same state but rotated along the incident direction over 
an angle 6, the elastic scattering angle. For most of the final states \n>, the transverse 
momentum of produced particles, p?, relative to the incident direction is sharply cut off, 
and its average increases slowly with s: the distribution in rapidity-pr space resembles 
that of a uniformly filled cylinder, sometimes called a "Wilson-Feynman liquid", with 
short-range correlations only between the hadrons. For such a configuration, it is easily 
verified that the inelastic overlap function, and thus Im^4 e /, will fall-off as an exponential 
in t, at small \t\, with a slope determined by the mean number of particles produced 
and by their (p\). For example, in a model where particles are produced independently, 
one finds |k| (see also fl65| ) 

B(s) > constant + {n)/(p 2 T ). (37) 

Consequently, B(s) grows with energy like (n), the mean multiplicity of produced 
hadrons if (p^) is constant. This explains the shrinkage of da e i/dt. For this estimate 
phases of the multiparticle amplitudes are neglected. The phase of the amplitude is 
related to the position in space-time where the particle is produced [66], and is unknown. 

Writing (n) = u Ay = u ln(s/s ), we see that B(s) depends on the particle density 
in rapidity space in inelastic collisions and by the variance of the transverse momentum 
distributions. In more rigorous calculations, the second-order transverse momentum 
transfer correlation function^ enters in equation (|37D instead of (pfi) ||4~8l, |6"7| . 

This result is generic and valid in a wide class of models (see e.g. |68| ). In processes 
where (p^) is larger than a soft scale, or large compared to \t\, the second term on the 
right-hand side of equation (|37| ) will be unimportant and shrinkage will either be small 
or absent. This most likely happens in (quasi-)elastic processes where a large scale can 
be identified ("hard diffraction"). 

Whereas the overlap of the amplitude of two "Feynman- Wilson liquids" will be 
negligible at large t, one realizes easily that hard jet emission will contribute to non- 
zero values of the overlap function at large t. This is the basic reason for the importance 
of very large t scattering and its connection with perturbative QCD. 

In a general collision process, and 7*p in particular, both u and (pf.) can be 
expected to be process- and (perhaps) energy- dependent. There is no sound reason to 

% In pQCD "ladder-language", this is the correlation between neighbouring propagator transverse 
momenta. 
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believe that these quantities, and thus the intercept and slope of the dominant Regge 
trajectory, are universal. 

The arguments given show clearly the connection between properties of the final 
states and Regge trajectory parameters for diffractive scattering. Since the relevant 
dynamical quantities, here uo and (p?), are clearly identified, generalization beyond the 
Regge framework is, at least conceptually, simple to understand. 



5.2. Inelastic diffraction as a regeneration process 

The possibility of inelastic diffraction has been predicted in the seminal papers by 



Feinberg and by Good and Walker [41]. Consider a projectile (hadron, real or virtual 
photon, etc.) hitting a target at rest. The projectile, being composite, can be described 
as a quantum-mechanical superposition of states containing various numbers, types 
and configurations of constituents. The various states in this superposition are likely 
to be absorbed in different amounts by the target. As a result, the superposition of 
states after the scattering is not simply proportional to the incident one. Hence, the 
process will, besides elastic scattering, also lead to production of inelastic states with 
the same internal quantum numbers as the projectile. This is the fundamental basis for 
inelastic diffraction and requires little more than the superposition principle of quantum 
mechanics, unitarity and the coherence condition, equation (p4). 



Assume that the projectile, \B>, at a fixed impact parameter (b) from the target 
is a linear combination of states which are eigenstates of diffraction 

\B) =J2C k \V k ), (38) 

k 

ImT|%) (39) 

where ImT is the imaginary part of the scattering operator and the (real) eigenvalue 
tk is the probability for the state \^k) to interact with the target. The eigenvalues 
or absorption coefficients tf. of course vary with b. The states are normalized so that 
< B\B >= J2k \Ck\ 2 = 1- From equation (|38|) and equation ( |3~9|) one easily derives 

da T /d 2 b = 2(t), (40) 

da el /d 2 b = (t) 2 . (41) 

The cross section for inelastic diffractive production, with elastic scattering removed, is 
^Sff /d 2 & = (t 2 ) - (t) 2 . (42) 

The brackets < ■ • • > denote an average of tk or t 2 ., weighted according to their 
probability of occurrence, \Ck\ 2 , in \B). We note the important result that inelastic 
diffraction is proportional to the variance in cross sections of the diagonal channels. 
Elastic scattering, on the other hand, is proportional to their mean value. Equations 
further imply the upper (Pumplin) bound |6£ 



°m(P)+or el (b)< l -a tot (b). (43) 
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From equation (|42| ) follows that, if the variance is zero (e.g. when all states are 
absorbed with the same strength) there is no inelastic diffraction. Diffraction will be 
strongest in regions of 6-space where absorption shows the strongest variation i.e. at the 
"edges" of the target: hadronic inelastic diffraction is more peripheral than the elastic 
process which is largest at small b. Further, in the case of complete absorption at a 
given b, inelastic diffraction vanishes at the same b. 

Note that for virtual photon scattering, the purely elastic reaction can be neglected. 
In this case, the term < t > 2 in equation ([42]) is absent. For real and virtual photon- 
hadron interactions, very little is known experimentally about the impact parameter 
profile. It requires a measurement of the t-dependence over a wide range in t. For 
elastic p production it was studied for the first time in [70|, following the method of 
Amaldi and Schubert f34j] . 

As remarked in ]7l|], "The increase of <tt (in DIS) with energy occurs because some 
regions of impact parameter are changing from grey to black and regions at larger b are 
going from white (no absorption) to grey. However, the region where absorption shows 
the strongest variation, and which contributes to diffraction, grows less rapidly than 
those 6-regions giving elastic and highly inelastic scattering. This would explain the 
observation that the inclusive diffractive cross section grows less rapidly than expected 
from Regge arguments (cfr. figure |T7|). Regge theory is indeed expected to hold for those 
regions in b where the absorption is weak. Regions of large absorption then correspond 
to multiple Pomeron exchanges." 



6. A generic picture of high energy collisions 

Well before the advent of QCD, and inspired by the ideas of Ioffe, Feinberg, Gribov, 
Pomeranchuck and others, a basic, although semi- quantitative, understanding of the 
space-time evolution of a high-energy scattering process was developed [^(].[7^.[73|]. It 
testifies to the profoundness of these ideas that, in spite of major developments in the 
field of strong interactions, the physical picture then developed still remains valid to a 
very large extent. Perturbative QCD has allowed us to clarify many issues, and produce 
crisp quantitative predictions in some cases, but the ideas then formulated, reaching 
well beyond pQCD, continue to be of great value. They provide a view of the collision 
dynamics which is simple enough to help develop intuition, provide physical insight and, 
hopefully, inspire new directions for future experimental research. 



6.1. Ioffe time 

It was first observed in QED |74| that photon emission from electrons propagating 
through a medium occurs over distances which increase with energy. In their 
seminal paper [75| Gribov, Ioffe and Pomeranchuk demonstrated that at high energies 
large longitudinal distances, now usually referred to as coherence-lengths, l c , become 
important for any kind of projectile, including virtual photons, when considered in the 
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rest frame of the target. The typical time involved is 0(E/n 2 ), with E the energy of 
the projectile and \i a hadronic mass scale, and basically follows from relativity. For 
DIS, Ioffe JFB| demonstrated that the longitudinal distances involved, measured in the 
target rest frame and in the Bjorken limit, are growing as 

lc oc — , (44) 

m,N x 

where is the target mass. It should be noted, however, that scaling violations, which 
are especially strong at small x, modify equation (|44]) and reduce the value of l c f7"7j] . 
In addition l c depends on the polarization of the virtual photon. 

The value of l c becomes large for small x or large W . At HERA, for Q 2 = 10 GeV 2 , 
the x values range between 1CT 2 and 10~ 4 and l c corresponds to distances of up to 1000 
fermi. Pictorially speaking this means that partonic fluctuations of the virtual photon, 
the Fock states, are long-lived and travel a substantial distance before interacting. 



< n > 



Figure 21. A high energy interaction in the parton model. From BS 



6.2. The Gribov-Feynman parton model 

The parton model views a high-energy interaction of any projectile, particle a, with 
a target, particle b at rest as follows (see figure ^TJ). The fast hadron fluctuates into 
point-like partons: quarks and gluons. The fluctuations have a lifetime t oc M% before 
interaction with the target occurs. During this time the partons are in a coherent 
state which can be described by means of a wave function. Each parton can, in turn, 
create its own parton cascade, each creating (n) partons, resulting eventually in the 
emission of a total of iV soft partons ( "wee" partons in Feynman's terminology |]72|| ). 
The latter should not be confused with the partons of pQCD. They are non-perturbative 
("dressed") objects due to the long time-evolution of the cascade and have acquired a 
large transverse extension. They interact with a target with a large hadron-like cross 
section. 

For a highly virtual photon, the cascade starts with a qq fluctuation (or dipole) 
of small transverse extension and is followed by an initial evolution stage where the 
strong coupling remains perturbative and calculable. However, the non-perturbative 
end of the cascade is likely to be similar to that originating from a high-energy hadron. 



Diffractive Scattering 



28 



Since gluons rather than quarks will be the dominant component of the cascade, and 
since gluons carry a larger color charge, it should not have been a surprise to find that 
the interaction, when viewed in an infinite-momentum frame, is driven by the gluon 
constituents in the proton. The same holds for "diffractive" structure functions of the 
proton or of the Pomeron. 

As argued by Gribov JT3J, a fast projectile can interact with the target only through 
its wee component. Indeed, the cross section of interaction of two point-like particles 
with large relative energy, y/s ab , is not larger than ttX 2 ~ l/s a b ~ exp(— r)ab) (A is the 
wavelength in the cm. frame of a, b, r]^ is the relative rapidity). Thus, only slow partons 
of the projectile are able to interact with a non-negligible cross section. Since there are 
iV wee partons in total, the interaction cross section is proportional to the probability 
that at least one wee out of N interacts with the target. For small N this is proportional 
to N. 

The interaction can also be viewed from the rest-frame of the projectile, or from any 
collinear rest-frame. The distribution of the wee partons in the rest-frame of particle 
b is, according to the above arguments, solely determined by particle a and does not 
depend on the properties of particle b. On the other hand, in the rest frame of particle 
a the distribution is determined by the properties of particle b. This is possible only if 
the distribution of partons with rapidities r] much smaller than the hadron's rapidity, rj p , 
does not depend on the quantum numbers and the mass of that hadron. It follows that 
the distribution of the wee partons with 77 -C rj p should be independent of the projectile 
and target, i.e. be "universal". Indeed, in the cascade the memory of the initial state is 
lost after a few steps only, if it resembles a Markov process. 

The fact that the wee-parton component of any hadron is independent of the 
hadron itself, explains semi-qualitatively why hadronic total cross sections show a 
"universal" energy-dependence at large s, as discussed in section |3j]. In addition, 
if the interaction between wee's is effectively short-range in rapidity (implying that 
the produced hadrons show short-range rapidity correlations), hadrons produced in 
regions of rapidity sufficiently far from target and projectile will also show "universal" 
properties. 

6.2.1. Shrinkage. Consider the interaction in the impact-parameter plane, figure 
In each step of the cascade the newly emitted parton acquires a certain amount of 
transverse momentum, kr- If the emission is purely random in /c^-space, the last parton 
in the cascade will, as the result of a random walk in impact-parameter space, have 
moved a distance b 2 N from the origin. On average, and for a completely random process 
{which precludes any kind of px- ordering of the emissions), one has 

( ^ )oc W) N = W) lns/s °' (45) 

for N cx In s. In this simplified picture, the (transverse) growth of the interaction region 
with energy is thus the result of a diffusion process. It is represented as the shaded area 
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in Fig. Thus, u /(k^) can be identified with a' F in equation (|T4|). The argument 



based on the overlap function, section [5.1|, leads to the same result but is more general. 

It is superfluous to mention that wee-parton properties, and their interactions, 
cannot be calculated in pQCD, neither for hadronic collisions nor for deep-inelastic ep 
scattering, since they are associated with long-wavelength fluctuations of the color fields. 
For DIS, this ignorance is parameterized in the parton distributions at a small scale. 
However, in the small-x region, the wee partons are equally important in both types of 
scattering processes. 



6.2.2. Rise of the total cross section and a F (ff). Consider figure [21]. Since each parton 
in the parton cascade can form its own chain of partons, and so on, this multiplication 
process will generically (but not in detail) lead to a total mean N ~ e^ n \ if (n) is mean 
multiplicity in a single chain ||68|| . With (n) oc ujq Ins one finds: 

iVocs" . (46) 

This can be rewritten in a frame-independent form 

1 f sY° 

o~t = cr (projectile) x o~q (target) x — x I — . (47) 

s V s o/ 

The "impact-factors", do, are particle-specific but independent of s. Equation ([47]) 
"explains" the power-law (or Regge) behaviour of o"t- 

For a collision of a small-size (in 6-space) virtual photon with a proton, Q 2 larger 
than a few GeV 2 , the evolution of ujq with Q 2 is calculable in pQCD. This is one of 
the major theoretical advantages of deep-inelastic scattering over soft hadron-hadron 
interactions. Evidently, in DIS, the role of s is taken over by W 2 or 1/x. 

From equation fl4"T| ) we see again that a P (0) — 1 in Regge theory has to be interpreted 
as the wee-parton density in a parton cascade. This result is generic. However, since the 
detailed process-specific dynamics of the parton cascade (DGLAP [ff8| , BFKL |79fl ,. . . ) 



will influence the evolution of o>o, we may conclude that a "universal" Pomeron trajectory 
with process-independent parameters does not exist. 

The power-law form, equation ((47]), is a result typical for a self-similar (fractal) 
branching process with fixed coupling constant and u>o is related to the fractal dimension. 
Early pre-QCD examples can be found in [B0|j. For a running coupling constant, the s 



dependence is generally less strong, but faster than any power of In s. 

6.2.3. Total cross sections, diffraction and wee-parton multiplicity. Suppose the 
projectile is a superposition of states with, at given impact parameter b, n wee partons, 
each of which can interact with the target with a probability f(b). If the structure of 
the target is ignored, we have (for brevity, we omit the argument b in the following) 

a tot = J2 a tot(n)P(n); (48) 

where P{n) is the probability that the cascade has produced n such partons. Using 
conservation of probability (or unitarity) we find 

a tot {n) = 2T el {n); T el {n) = I - - a in {n); a in (n) = 1 - (1 - /) n . (49) 
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The last equation in (pE9| ) is the probability that at least one out of n partons interacts 
with the target. 

The previous equations can be compactly expressed in terms of the generating 
function of P(n), E(z) 



= P(n)(l + z) n , (50) 

a tot = 2 p (n) [1 - (1 - f) n/2 } =2-2 E(y/T=f - 1), (51) 

cT dlff+el =J2 P (n)[l-(l-f) n/2 } 2 , 

= 1-2S( > /W-1)+E(-/). (52) 

For the ratio of total diffractive (sum of inelastic and elastic) cross section to the 
total cross section, R(b), at fixed impact parameter, we obtain 

{h = <y_^f I± e L 1 1 -^LB r . (53) 

1 ' o-tot 2 l-S(VW-l) 1 ' 

In the case of total absorption, / — > 1, the ratio converges towards the black-disk limit 
of |, as it shouldF]. 

Assuming, as an example, P(n) to be Poissonian we obtain 

o-tot ~ ^/ < n >, (54) 

o-diff+ei « — < n 2 >= — [< n > 2 + < n >] ; (55) 

provided / or (n) or both are small enough. These conditions mean that multiple 
interactions with the target can be neglected, or that the partonic system hitting the 
target is sufficiently dilute and no saturation takes place. 

Equation (|>|) suggests a relation between the total cross section (or F 2 in DIS) 
and the mean parton multiplicity which was first tested experimentally in |H] and is 
illustrated in Fig. [22|. Using a Modified Leading Log (MLLA) pQCD expression for 
the energy-dependence of (n) (Eq. (7.32) in [0), an excellent description of the in- 
dependence of F 2 {x, Q 2 ) data at low x was achieved with two free parameters only. 

If saturation (parton-recombination) effects in the parton cascade in DIS happened 
to occur, we can expect that the similarity of the energy-dependence between mean 
particle multiplicities in e + e~ and F 2 will break down for very high multiplicity events. 
Given the present interest in this topic ||29||, a dedicated measurement of the In- 



dependence of semi-inclusive structure functions F 2 (x, Q 2 ), at fixed large final-state 
multiplicity n, and of its diffractive counterpart, might therefore be of considerable 
importance. 

Assuming a s e dependence of (n), we further see that ([54]) and (|55f) predict the 
energy-dependences o tot oc s e and (Tdiff+ei oc s 2t , the same as obtained in Regge theory, 

+ For DIS at very large W, it follows, quite remarkably, that for scattering on a very large fully 
absorbing nucleus, 50% of the total cross section will be diffractive, even when Q 2 is very large (but 
Q 2 /W 2 < 1). 
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Figure 22. Comparison of e + e data on average charged particle multiplicities 
versus yfs and the HERA low- a; data versus 2Qi/y/x, with Qi = 270 MeV, for 
Q 2 = 22 GeV 2 (ZEUS) and 25 GeV 2 (HI). The e+er multiplicity data (solid lines) 
are represented by curves resulting from a phcnomcnological fit of a second-order 
polynomial in Ins to e + e~ data. They are normalized to the F2 data for each Q 2 
bin separately. From pM. 



and thus show the same unitarity violating defects as mentioned in section |2.2.3| . To 
obtain a constant ratio, R(b), at each impact parameter, it seems unavoidable to include 
in the calculation the full multiple-scattering terms and possibly (so far unknown) parton 
correlations. 

The role parton correlations can be illustrated using the factorial cumulant 



expansion of H(z) (see e.g. [|3|, 



The cumulants K q are a measure of the correlations and identically zero for q > 1 if the 
partons are uncorrelated. The inelastic cross section can now be written as 



1 - »(-/) = 1 - exp \ -N f + J2 (-/)VffO K* ■ (56) 



Comparing equation ( |3"2"D with equation flSBP we see that the eikonal function Q(b) 
can be expressed in terms of the cumulant generating function In S(— /). This shows that 
not only multiple scattering contributions, but also parton-parton correlations (provided 
that K q 7^ for q > 1) contribute to the total and diffractive cross sections. Such 
correlations have not been explicitly taken into account, as far as we know, in present 



pQCD calculations of DDIS, with the exception of |85| using the concept of (Mueller) 
dipoles in onium-onium scattering. 
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6.3. Models for diffraction 

6.3.1. Diffraction and the parton model: the Miettinen and Pumplin paper. The first 
detailed calculations of hadronic diffraction in the framework of the parton model were 
presented in [j36|. This work, although 22 years old, remains of great interest and we 



summarize its main conclusions. 

It is assumed that the diagonal states (section |5?2|) |^ > are the states of the parton 
model, composed of quarks and gluons and a radiation cloud of wee partons. These 
states are characterized by a definite number N of partons with impact parameters 
&jv an d longitudinal momentum fractions, or rapidities, yx, . . . , ypf. 

Since there are parton states which are rich in wee partons, and others with a few 
or no wees, these states will interact with a target with very different cross sections. 
Hence, inelastic diffraction will be generated by the mechanism of Good and Walker. 
The fluctuations in the interaction probabilities tk (equation (^)) arise from fluctuations 
in the number of wee partons, fluctuations in yi and from fluctuations in bi. 

Assuming uncorrelated wee partons, and fitting all free parameters of the model 
to o c \{pp) and a^ipp) at ^fs = 53 GeV, the calculated total inelastic diffractive cross 
section was found to be in very good agreement with data. The yi fluctuations contribute 
little (about 10%), whereas fluctuations in bi and in parton number each account for 
about 45% of crP D . Also the forward value and the slope of the t distribution are 
correctly predicted. This is a non-trivial result since the calculated (and measured) 
slope B pa 6.9 GeV~ 2 is only about half that of elastic scattering B pa 12 GeV -2 . 



Interestingly, as seen from figure [23], the small \t\ dissociation is dominated by the 



large and very steep (slope ~ 12.2 GeV 2 ) contribution due to the parton-number 



fluctuations, see also equation (|55f ). The bi fluctuations, on the other hand, dominate 
at large \t\. 

The ZEUS collaboration recently presented new measurements, shown in figure |24| 
(see [0), of the t-slope in diffractive DIS, using their Leading Proton Spectrometer 
(LPS). The slope has a value B = 6.8 ± 0.6 (stat)+^; 2 (syst) for 4 < Q 2 < 100 GeV 2 , 
Mx > 2 GeV, £ < 0.03. Some evidence for shrinkage is seen but no dependence on Q 2 . 
The value of the t slope is strikingly similar to that in the pp data. 

6.3.2. Modern QCD models of diffraction It is evident that the Miettinen-Pumplin 
(MP) model grasps the essential physics which remains valid in the context of DDIS at 
HERA. In |^(| ad-hoc assumptions were needed to build a model of the hadron Fock 



states. In DIS the light-cone "wave functions" of the lowest-order 7* Fock states (qq, 
qq + gluon) are known |5(||88[] and quantitative results can be obtained. Nevertheless, 



the interaction of the wee partons needs to be parameterized empirically, as it must be 
for soft hadron-hadron collisions. 

The presently popular models for diffraction in DIS have been reviewed in 0. They 
use the same basic concepts discussed in previous sections under various disguises. The 
most successful of these are merely modernized versions of the Aligned Jet Model |Kj 
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Figure 23. Differential t distribution, da dt , and model calculations. The 
decomposition is shown of the full cross section into contributions due to fluctuations in 
the number (N), rapidities (yi) and relative impact parameters (pi) of the wee partons. 
The iV-fluctuation dominates near t = 0, and the bi fluctuation component at large t. 
From [ p6[ . 
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Figure 24. ZEUS: the t distribution measured in D1S, using the Leading Proton 
Spectrometer [p7f . The slope has a value B = 6.8 ± 0.6 (stat)lj'y (syst) 
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which pre-dates QCD. 

Considered in the target rest frame, the Fock-state wavefunction of the partonic 7* 
fluctuation carries the information on the virtuality of the photon and further depends 
on its transverse size and the fractional momenta and masses of the partons. In the 
simplest case of a qq fluctuation, or qq dipole, the wave function is then convoluted with 
the amplitude for the elastic interaction of the colour dipole and the target hadron. 
At t = 0, this amplitude is determined by the cross section for the scattering of the 
dipole with the target, a(g). It is assumed to be independent of Q 2 , in accord with 
the Gribov-Feynman argument of wee-parton scattering and short-range order in the 
cascade, but depends on x. 

Consider, as an example, the very successful saturation model of Golec-Biernat and 
Wusthoff (GBW) |H| which expands on much earlier work ||88|| . The physical picture 
is that in which, in the nucleon rest frame, a photon with virtuality Q 2 , emitted by a 
lepton, dissociates into a qq pair far upstream of the nucleon. This is then followed by 
the scattering of the colour dipole on the nucleon. In this picture, as also assumed in 
the MP model, the relative transverse separation q of the qq pair and the longitudinal 
momentum fraction z of the quark remain essentially unchanged during the collision. 



The 7*p cross sections take the following factorized form f88|,|89| 

(T T Ax,Q 2 ) = / d ^/ dz\^ T>L (z,g,Q 2 )\ 2 a(x,g), (57) 

where *&t,l is photon wave function of transversely (T) and longitudinally polarized (L) 
photons. 

In (0), all Q 2 dependence is contained in the Fock-state wavefunction, which 
further depends on the flavour and mass of the partons. The W- or x-dependence of 
&t,l{x, Q 2 ) is solely determined by that of a (x, g). The latter is the principal quantity 
in the s-channel description of diffractive scattering. Once the dipole cross section is 
known, (|5"7j ) enables a parameter-free calculation of the proton structure function at 
small x. In our simple picture, we may interpret it as an effective cross section, the 
product of the wee-parton flux with the single wee-parton nucleon cross section. 

Although the impact-parameter dependence of a (2, g) is not explicitly considered 
(only its average enters in (|57D ), this is clearly of great interest for the t-dependence of 
the diffraction [|7(]], and needs to be studied further. 

Turning to diffraction, the differential cross section at t = takes the form 



dt 



16 

i=0 



^ Jd 2 g^dz\^T, L (z,Q)\ 2 (T 2 (x,g). 



(58) 



The form of (|58| ) differs only from ( |57D by the substitution a(x, g) — > <r 2 (x, g), in accord 
with the general formula (^) 

* In the dipole formulation of DIS, and contrary to hadron diffraction, even inelastic DDIS is considered 
to be purely elastic: the dipole states (qq) and higher-order Fock-states qq + gluons are assumed to be 
orthogonal eigenstates of the diffractive T-matrix, and no regeneration (mixing of the states) occurs. 
If these states are not orthogonal, they will regenerate and thus add an additional contribution to the 
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Comparing to the MP-model, we see that the relative impact parameter and 
rapidity fluctuations are included here through the photon wave function. The 
important parton-number fluctuations, which also depend on parton-parton correlations, 
however, are not explicitly considered. 

The energy-dependence of a (x, g) follows from the fact that, in low-a; DIS, the 
perturbative evolution of the qq dipole results in further "hard" parton multiplication 
which increases also the wee-parton flux and thus the total cross section. Indeed, due to 
the bremsstrahlung nature of soft gluon spectrum oc dz g / z g (z g is the momentum fraction 
of the photon carried by a gluon) Fock states with n such gluons give a contribution 
oc In (l/x) n to the total photoabsorption cross section, which can be reabsorbed into an 
energy-dependent dipole cross section f88 |. For example, in the DGLAP aproximation, 
summing over all n produces the well-known exp [2^/ln (1/a;) In (l/a s (Q 2 ))] increase of 
the 7*p cross section and "standard" scaling violations. 




Figure 25. The dipole cross section a(g) for various values of Bjorken-x. The GBW 
parameterisation, Eqs. 60 and O, with A = 0.29 is shown as the solid curve. The 



dashed lines correspond to the 1/x dependence given by Eq. |62j with parameters n/ = 3, 
K = 0.288, A = 1.03 GeV, taken from (9fJ. The dot-dashed lines show 0.05 x <r 2 (g). 



The qq dipole-proton cross section a(x, g) has to be modelled although it is known 
in the perturbative limit of very small dipoles and related there to the inclusive gluon 

diffractive cross section, presently neglected. 
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distribution xg(x,fi 2 ) of the target |9~T| . 



7T 



(59) 



valid at small g; C is a scale parameter. For a (ggg)-dipole system, a(x, g) is roughly 
a factor 9/4 larger. This explains the predominant role of gluons in low-x DIS. In the 
GBW model, the effective dipole cross section is taken to be of the form 



a(x, g) 



00 



exp 



Q 



AR 2 {x) 



where the x-dependent radius Ro is parameterized as 
1 . /^„\A 



Ql p 

V x 



(60) 



(61) 



with Qo = 1 GeV. The parameters o"o = 23 mb, xo = 3 . 1CT 4 and A = 0.29 have been 
determined by a fit to data on Fi |29j . As seen in figure |2^, the dipole cross section 
saturates at a value o"o for large-size dipoles where it is entirely non-perturbative. Also, 
as x — ► 0, saturation sets in at decreasingly small transverse sizes, and the contribution 
from large-size dipoles becomes more important. 

Since equation fl58l) depends on the square of a(x,g), it follows that still larger 
sizes are involved in diffraction than those dominating the total cross section: non- 
perturbative soft physics is of even greater importance in DDIS (see the dot-dashed 
lines in Fig. |25[). Saturation effects are therefore predicted to be more important than 
in inclusive DIS. Because of the saturation property of (0), nearly the same dependence 
on x and Q 2 of DDIS and DIS is found, thus giving a natural explanation of the constancy 



of their ratio as mentioned in section fO 



In the GBW model, two essential scales appear: the characteristic transverse size 
of the qq dipole oc 1/Q, solely determined by the 7* wave function, and Rq(x). Naively, 
1/Rq(x) can be interpreted as the mean number of soft partons in the cascade; Rq(x) 
is their mean relative transverse distance and Q Rq(x) = 1 defines a critical line. For 
1/Q <C Ro(x) the partonic system is dilute, for 1/Q ^> Rq(x) the system is densely 
packed and multiple scattering and parton-interactions become important. 

It is interesting to note here that the fitted value of A (~ 0.29) is quite close to 
that derived from the cm. energy-dependence (v^) of the mean particle multiplicity 
in e + e~ annihilation, where it is found that < n >~ s 25 provides a reasonable fit of 



the data [92]. Remembering the striking analogy discussed in section |6.2.3 , we have 



also plotted in figure |25| expression (|60"D wherein 1/ R^(x) in Eq. (|61| ) is replaced by that 
of the mean soft gluon multiplicity in a gluon jet with energy-squared oc 1/x, as given 
in 



90 



KQ 2 Q y- a ^e XV [2C^ + 5 G 



with K an overall normalization constant, C 



Sc(y) = 71 1 



2a 2 C 2 + §M2») + 2] 



C 2 
+ — 

y 



= ^/4N c /(3 , and 
azC 2 - a -^[M2y) + l] 



(62) 



(63) 




Figure 26. Anomalous multiplicity dimension to different orders in The chosen 
example is Nf = 4 and A = 0.22 GeV. Anomalous dimension N' g /N g as a = 3a s /2w 
(solid line), \fK{l— 2a\y/a) (crossed) and yJ~E(l — law/Ik— 4ci2<5) for the parton equation 
(dashed), dipole equation (dotted) and generalized dipole equation (dash-dotted). The 
constants are a\ = 0.297, ci2 — —0.339. For a detailed discussion see |94| ]. 



Here e y = y/l/x/A and further (3 = (11N C - 2n f )/3, pi = [17iV c 2 - n f (5N c + 3C F )]/3 
N c = 3 is the number of colours and Cp = 4/3. The numbers are tabulated in 
A is the QCD scale parameter and rif the number of active flavours. 



The dashed curves in figure |25| show the dipole cross section as obtained from 
equation (|62|) . It is essentially indistinguishable from the GBW parameterisation for 
x = ICT 2 — 10~ 3 , but differences become noticable at smaller x. This follows from 
the fact that, due to the running of a s , the multiplicity grows slower than a power in 
1/x and "saturation" is delayed in comparison with Eq. (|61]), the latter being a result 
characteristic of a cascade process with a fixed coupling constant. 

The results shown in Fig. ^5] imply that the ansatz in Eq. fl62|) will lead to an 
equally satisfactory description of F 2 and as was obtained in the original GBW work. 
However, the parameterization (0), contrary to (0), involves no free parameters, apart 
from Q and the normalization constant K which was taken from a fit to e + e~ data PC 



In particular, the important parameter A follows in the former case from theory. 

If 1/Rq(x) is interpreted as the mean number of soft partons confined in the target 
within a transverse surface of radius R, it is evident that the GBW parameter A can be 
identified with the anomalous (multiplicity) dimension, 7, of the parton cascade (see 
for a recent review and further references) which is calculable in pQCD. 

With L = ln(s/A 2 ) (s is the relevant energy-scale squared, 1/x for DIS), 7 is in 
general defined as 

ry 

N g oc exp / 7(2/) dy' . (64) 
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where N g is the mean gluon multiplicity. In pQCD, for a time-like cascade, it is equal 
to the logarithmic derivative of iVg with respect to L and given by |90 | 

7 = = y/E(l- 2a x V6c-Aa 2 a + 0(a 3/2 )) (65) 



The first term, \/<f = a/ N c a s /2ir, is the leading-order term. However, 7 is non-linear in 
a and decreases with increasing s. Figure [26| shows a plot of the anomalous dimension, 
calculated to different orders in a and for various cascade schemes (parton-cascade and 
Lund-dipole pictures) as discussed in f94|. It is seen that 7 decreases with s due to 
the running of a s . Consequently, N g increases slower than a power of s. A power-law 
dependence is obtained if a s is kept constant. Taking, as an example, a s = 0.2 in ( |65f) 
yields A = 0.30 at lowest order. 

The relation between the GBW parameter A and the multiplicity anomalous 
dimension has also been derived in the framework of the Balitsky and Kovchegov 
modified BFKL equation (see |96| for details and references) with the result A ~ 6a s /7r. 
For a s = 0.2 this gives A = 0.38. Note, however, that ( |65|) is a polynomial in ^foT Sl 
whereas the previous expression is linear in a s . 

An further important result of the GBW-model is that the diffractive structure 
function is found to obey a Regge-like factorization property (except for /3 — > 1 where 
higher-twist contribution from longitudinal photons dominate) with the dependence 

oc (1/£) 1+A . This corresponds to an effective Pomeron intercept a P (0) ~ l+A/2 f97|j . 
This is a highly revealing result, demonstrating, on this specific model-example, the 
generic property that the growth of the cross sections with energy and the proton 
energy-loss spectrum (or "Pomeron flux") are closely connected (cfr. section |4.3.1| ) 
and determined by the multiplicity anomalous dimension. 

The above considerations suggest an extremely simple picture (see also [^] p. 8) 
for the x and Q 2 dependence of F2 and _F 2 D . At low x, the target is populated with 
a number of partons proportional to N g confined within a transverse area 7r R 2 . Since 
the area "scanned" by the virtual photon qq dipole is proportional to 1/Q 2 , the number 
of partons with which it can interact is proportional to - = N g ^2 °^ ^ ^2 ■ One can 
therefore expect that the total cross section will depend only on r and not on Q 2 and Rq 
separately, i.e. exhibit "geometrical scaling" J99|. This follows already from dimensional 
arguments but also agrees with the "universality" hypothesis, advanced in |100|] , that 
the physics should depend only on the number of partons per unit of rapidity and per 
unit of transverse area. If 1/r is sufficiently small, multiple scattering effects can be 
neglected. 

To see the influence of multiple scattering, we return to formulae (|5lT) -(|52"[) which we 
now apply to the cross section of a dipole of fixed size g interacting with the target. We 
further assume that parton correlations can be neglected. In that case, the generating 
function E(z) is that of a Poisson distribution E(z) = exp (< n > z). Provided that / 
is small enough we obtain 

a dipoie ~ 2 (l-e-? <n>f ) (66) 
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ag55 - Kl-2c-*<->' + c-<->' (67) 

Equation (|66| ) is precisely of the GBW (eikonal) form (|60| ) if, following the previous 
arguments, < n > f is identified with 1/r, the effective number of partons "seen" by 
the qq dipole. 

The above formulae invite further comments on the meaning of the term 
"saturation". The form of equation (|5SD follows from that of the generating function 
which includes the full Glauber-Mueller multiple scattering series which leads to a 
levelling off of the dipole cross section. Only for a very dilute parton system, or for 
a very small dipole can these additional terms be neglected. On the other hand, parton 
recombination effects, when they occur, will induce a weaker 1/x dependence of < n >, 
compared to that given e.g. by Eq. (|62D. Although our simple semi-classical picture 
therefore suggests two distinct origins of saturation, it is not clear if such a distinction 
is physically justified in more rigorous treatments of the dynamics. 

The model results, discussed previously in the target rest frame, can be translated, 
at least in leading twist, in terms of diffractive parton densities in an infinite momentum 
frame. The diffractive structure functions are then expressed as the convolution of 



"diffractive" parton densities for the proton with parton cross sections ||101]| . The 
evolution with Q 2 at fixed x is the same as that of F 2 (x,Q 2 ). In DDIS, these scaling 
violations affects the (3 (or Mx) dependence of the cross section but not the dependence 
on £ ||102|| . However, unlike the case of fully inclusive cross sections, the diffractive 



structure functions are no longer universal. In particular, they cannot be used directly 



for hadronic interactions 97 . 



7. Summary and outlook 

Over the last decade, the subject of diffraction has become one of the very active fields 
of experimental and theoretical research in QCD. The revival is, by large, due to the 
extremely varied experimental programme made possible at HERA and at the Tevatron. 

In this paper, we have attempted to describe, mainly in qualitative terms, the close 
relation between the dynamics of total cross sections and diffraction in hadron-hadron 
collisions and in deep-inelastic 7*p scattering. This inter-relationship is ultimately 
a consequence of the fact that the bulk of the total and diffractive cross section is 
dominated at very high energy by the wee components of the target and projectile's 
wavefunctions such that long-distance physics plays a very important role in both. 

We have argued that the physics can be understood on the basis of a surprisingly 
small number of dynamical ingredients such as the anomalous multiplicity dimension of 
parton cascades, 7, which not only determines the rise with energy of the cross sections 
but also the spectrum of the elastically scattered proton in DDIS. 

Our discussion of the overlap function illustrates that the small \t\ behaviour of 
the quasi-elastic processes is also determined by 7 and by the transverse-momentum 
transfer correlation function. These ingredients suffice for a basic understanding of the 
degree of "shrinkage" of the forward diffractive peaks in soft as well as in hard processes. 
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The remarkable recent theoretical progress in DDIS is a consequence of the 
fortunate circumstance that perturbative QCD is able to make reliable predictions for 
the partonic fluctuations (Fock states) of a virtual photon, and for the subsequent 
development of these states into a parton shower or radiation cloud, at least in the 
earliest perturbative phase of the evolution. For a strongly bound system of large size, 
e.g. a hadron, which is less well understood, such perturbative techniques are not 
available. 

Much of the present phenomenology of diffraction can be understood from the 
properties of the 7* Fock states. 

• The Mx distribution for the lowest-order qq dipole state, and transversely polarized 
7*, has the form l/[m 2 (Q 2 + M x ) 2 } (m,f is the quark mass). Extra soft gluon 
emission, with a spectrum dz g /z g directly leads to the much weaker M x dependence 
dM x /(M x + Q 2 ) in the so-called triple- Pomeron region [Q. Since the invariant 



mass of the diffractive system will remain almost unchanged for small-t scattering, 
this is also the distribution of the experimentally measured Mx- A similar argument 
was used in |58j to explain the \jM\ dependence of hadronic diffraction at large 
Mx- The gradual transition from a steep \jM x to a \jM\ dependence is a 
consequence of a change in the mixture of Fock states as Q 2 and/or W change. 

The average transverse size of 7* fluctuations relevant for elastic vector meson 
production, the so-called scanning radius PD[ , is estimated to be ~ C /{my + Q 2 ), 
with C ~ 2 (C ~ 6) for longitudinally (transversely) polarized 7* |[70| . This 
follows almost directly from the form of |^t,l| and from that of the vector-meson 



wavefunction. The elastic vector meson data (see e.g. figure |10] in section |3.3.2j ) 
show that scale Q e ff = Q 2 + my is indeed the dynamically relevant observable. For 
further discussion on this point we refer to |53] . 



• For high-mass diffraction, the partonic fluctuation of the 7* has the colour-topology 
of a gluon-gluon dipole in a colour-single state ||88|| . High- mass diffraction therefore 
opens the possibility, not yet fully exploited, to study the fragmentation of colour- 
octet sources, in much the same way as with qq + gluon three-jet events in e + e~ 
annihilation. 

The development of the radiation cloud initiated by the virtual photon is a cascade 
process, whereby the virtuality of the system is gradually degraded and a system of 
"perturbative" partons created. However, once the virtuality of the partons has reached 
values for which the strong coupling is no longer small, the cascade will continue into 
a non-perturbative region which is not under theoretical control. This corresponds to 
a regime in which the non-perturbative off-springs have acquired transverse dimensions 
comparable to the size of the target proton. It can be assumed that they will interact 
with the target as dressed objects, with a large hadron-like cross section. This and the 
variation in absorption is, by the mechanism of Good and Walker, the cause of shadow- 
scattering and diffraction, not only in DIS but also in hadron-hadron scattering. 

At large energy, the end of the parton cascade will show "universal" properties, 
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independent of the parton system which initiated the cascade. This, in turn, leads 
to expect factorization of the type often assumed in Regge theory and also found 
experimentally 

The simple picture described here suggests further experimental work in different 
directions. We end by listing only a few examples. 

• Studies of leading-proton and leading-neutron production in DIS and photoproduc- 
tion, also outside the diffractive region, combined with the many existing hadron- 
hadron data, should allow to test Regge-type factorization or provide evidence for 
factorization breaking. The latter is expected at low to moderate values of Q 2 . Ev- 
idence for possible long-range correlations between leading baryons and "central" 
hadronic activity (multiplicity, transverse energy density, jets) should be searched 
for. 

• The running of a s and 7, see Fig. ^6], suggest to measure in detail the £ dependence 
of ~ (1/£) 1+A as function of W and Q 2 . Whereas the kinematical range of 
the HERA experiments may be too limited to reveal an expected flattening of the 
da/dMx spectrum, running-a s and possible parton saturation effects which affect 
the value of A may become visible at LHC energies. 

• Since parton saturation is most likely to occur for high parton densities, a 
dedicated measurement of the W^-dependence of semi-inclusive structure functions 
F^\x, Q 2 ), and of its diffractive counterpart, for large final-state multiplicity (n) 
events could of considerable interest. 

• The dipole cross section a(x, g) plays a fundamental role in many models. At 
HERA, it can be measured in elastic vector-meson production but needs precise 
measurements of the differential cross section over a large range in t. 
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